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The Kondo effect underpins a large body of recent developments in the physics of d- and /-band 
compounds. Although its charge analog is a rarity in solids, the recent observations of the charge 
Kondo effect and the consequent rise in superconducting T c encourage a search for other accessible 
systems. Motivated by the possibility of wilfully tuning the sign of the inter-electronic interaction 
in optical lattices, we study conditions for the elusive charge Kondo liquid (CKL) state to manifest. 
We propose that a combination of Feshbach resonances and sequentially controlled laser pulses may 
produce the CKL. We show that the CKL is never a stable ground state, appearing only when the 
ordered ground states are destabilized. Finally, we discuss interesting analogies with nuclear matter. 

PACS numbers: 03.75.Ss, 71.10.Fd, 71.10.Hf 



The Kondo effect has historically played a monumental role in solid state physics where quasi-bound state 
(singlet) formation involving a local moment screened by the itinerant electron spins leads a renormalized Fermi 
liquid (FL) with huge enhancements of quasiparticle masses, whence the name "heavy fermion" . Subsequently, more 
exotic, quadrupolar [2| and orbital [3| versions were invoked in other contexts. The destruction of FL behavior at 
quantum phase transitions near magnetic order in some /-electron systems is thought to be linked to "unbinding" of 
this Kondo quasi-bound singlet state However, the intriguing possibility of observing a charge analogue of the 
Kondo effect, requiring an attractive, local interaction, has not been addressed carefully. In practice, materials such 
as Ba\- x {K,Pb) x BiOz are modelled by variants of U < Hubbard model [H, B S 0] • 

Recent experimental realization Q of the charge Kondo effect (CKE), predicted theoretically some time back [loj ]. 
and the consequent rise of superconducting T c in doped PbTe [11| and SnTe [9(] , has re- ignited interest in the U < 
models. Here, chemical doping (< 2%) generates a narrow impurity band in the semiconductor band gap, facilitating 
use of a U < lattice models in the intermediate-to-strong coupling limit (|J7|/Wj mp > 1). The interesting issues for 
theory are thus: What are the precise conditions under which a charge-Kondo Fermi liquid state (the CKL state) can 
appear? Is it ever a stable phase at zero temperature? Might a quantum phase transition (QPT), reminiscent of its 
well-studied spin analogue, occur, and if so, how? What instabilities might one generically expect to arise near such 
a QPT? 

Advances in artificially engineered fermionic/bosonic optical lattices allowing wilful manipulation of parameters 
open up yet another route to test model Hamiltonian predictions. Such systems may hold more promise for CKE, given 
the paucity of real materials exhibiting U < 0. In contrast to real materials, these systems are free of inhomogeneities. 
In this context, various interacting models have already been realized [13, [Hj]. The physical conditions for this lie 
well within the regime of experimental techniques available. In practice, different hyperfine states of same or different 
atomic species (acting, as it may, as different fermionic/bosonic species) can be trapped and controlled independently. 
A particularly attractive feature is their ready tunability: it is even possible to choose the sign of the interaction in such 
systems; a mixture of two fermionic species interacting via an attractive interaction is achieved by forcing a mixture of 
two such atomic spin states through a Feshbach resonance, whence a bound state appears in the two particle problem. 



A p-wave Feshbach resonance [LJ] can even create a tunable asymmetry in the interactions, allowing one to access 



more complex "hidden" ordered states. Further, the periodic potential of each species of atoms was independently 
tuned in an optical lattice [IH . This could facilitate observation of the band insulator (driven by staggered periodic 
lattice potential) to Mott insulator (due to on-site interaction, U) transition, and possible emergence of correlated 
metallic phases sandwiched between them. 

Optical lattices thus provide a unique tool to address interesting questions posed above. A suitable model for a 
two-component fermionic system (having an attractive two-body interaction), with temporally separated laser pulses 
simulating inter-site hybridization between the two spinless fermionic species (a, b) is described by the Hamiltonian 



H = } j t ab (a\bj + h.c.) - U ab ^ n ian lb + A >,(n, a - n ib ) + ^ £ a n ia 



(1) 
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The hoppings are t aa {tbb) for a — a(b — b) hopping and t a b for a — b hopping between nearest neighbors. A 
is the "charge transfer" term related to on-site fermionic energies. In typical optical lattice systems, the model 
parameters expressed in terms of the recoil energy are, hopping t — ^^■(■^■) 3 / 4 e~ 2 v / ( v '° < ,E and local correlation 

U ~ a s fci(Vb/-E'ii) 3/ ' 4 uM- Here, the recoil energy, Er = h 2 k\/2m is typically of order a few micro Kelvin, 

(fci and Vq are the wave- vector and intensity of the laser respectively). a s is the s-wave scattering length, and the 
fermionic bandwidth, W « Er. 

Typically, t/H ~ 1 kHz and U/% ~ 20 - 30 kHz 0. So U/W ~ 0(2 - 3) in three dimensions. We envisage 
comparable values for the hoppings and U a b in our case. This is precisely the intermediate to strong coupling regime 
of the U a b < "Hubbard" -like model that we study below. Focussing solely on possible Kondo-like physics, we 
use dynamical mean-field theory (DMFT) to study the model. Its reliability for understanding quasi-local collective 
Kondo coherence in correlated metals is already established [16(. Here, we show how a combination of experimental 
manipulations in optical lattices facilitates observation of the elusive CKL state. 

Eq. P is a generalized negative-U Hubbard model (relabeling a — > Cf,6 — >cjj with "spin" -dependent hoppings, 
and a "magnetic field" . Given the isomorphism between the U > Hubbard model and the U < Hubbard model 
in a magnetic field (seen by performing a particle- hole transformation), the possibility of the charge Kondo effect 
(CKE) involving pair-pseudospin-quenching by itinerance manifests itself (the pseudospin t" — c\ a a° b Cib). On the 
other hand, lack of pseudospin SU(2) invariance of H (lowered to Z 2 xU(l)) generically favors ordered ground states, 
either a charge-density wave (CDW) or a superfluid (SF). A CKL state can then appear as a ground state only upon 
quantum melting of these ordered states, or for model parameters such that Tckl > Tcdw.sf- 

To proceed, we "rotate" the a, b Fermions: writing f a = (ua + vb)/\/u 2 + v 2 , fb = (va — ub)/^/u 2 + v 2 with 
u = \Jt aa l(t aa + t b b), v = y/ t h b /{taa + Ub) on the parameter curve t ab = y/t aa tbb leads to 



<hj> 



where t — {t aa + hb), A = Aq — uv and e' a — v 2 e a , e' b — u 2 eb. This is the U < Falicov-Kimball model (FKM) with a 



local hybridization [17|. If A is tuned to zero, this reduces to the pure FKM. Since the FKM can be exactly solved in 



D = oo |18(, we first solve Hp km and then consider the effects of t a b ^ V^aahb, A^O. 

We observe that [riifb, Hfkm] = for each i, implying a local U(l) invariance of Hfkm- This is exactly the 
condition for degenerate, resonant charge fluctuations between local configurations with rif b = 0, 1 to occur. Elitzur's 
theorem implies that both local excitonic {{f} a fib)) and local pairing ((// a // b )) averages rigorously vanish. Choosing 
the / Q -band density-of-states (DOS) to be a lorentzian with half-width W; i.e, po{E) = W/n{E 2 + W 2 ), the / a -band 
propagator within DMFT is exactly calculable as 01 Gf a {u) = + ^W^+m ' where n fb = J2iifibfib) wi th 

N the number of lattice sites. Correspondingly, the /f,-Fermion propagator exhibits an infra-red singular, branch cut 
form, Gf b {u>) ~ 6{w)\u>\~^ 1 ~ ri ^ at low energy. Here, 77=tan _1 {U a b/W) /it is the so-called s-wave phase shift of the f a - 
band Fermions induced by the strong scattering term {U a b) in Hfkm- The /(,-Fermion self-energy, Hf b {uj) ~ |ll;| 1 ~' 7 , 
implying that the FL quasiparticle residue Z — 0, and the symmetry-unbroken metallic phase is a non-FL metal. The 
underlying cause for non-FL behavior is explicit in the impurity limit of Hfkm-, where recoil- less (infinitely heavy) 
/b-fermions scatter the / a -band fermions: it is the X-ray edge problem, where such an infra-red singularity is known to 
occur. This implies that the local pairing susceptibility (the analogue of the "excitonic" susceptibility for U > 0) also 
diverges in the infra-red, X/ Q ,/ 6 M = / dte 1 "* (fj a f] b {t); fibfia{0)) ^ M~ (2?) ~'' } : coupling of / a -band fermions to this 
singular pair susceptibility then destroys FL theory. Within DMFT, this is an incoherent pair liquid state, where pair 
correlations diverge in the infra-red without SC order. We dub this state the incoherent Bose metal (IBM). The IBM is 
characterized by resonant charge fluctuations as described above, and alluded to in earlier work (lol . |7|| . The mixing of 
the degenerate single-particle and pair states reflect the superposition of the "emission" and "absorption" initial states 
of the (presently, attractive-U) X-ray edge problem. Since the /{,-fermions are rigorously localized, the "Fermi surface" 
now wholly corresponds to the dispersive / a -band quasiparticles, with a smeared dispersion reflecting unusually broad 
pseudoparticles, even at T = 0. Depending on {t aa +tbb)/U a b and along the curve t a b — y/t aa tbb, we generically find a 
first order transition near half-filling, where the total fermion density, n, jumps from n < — {n/2) + ^tan~ 1 {U a b/2W) 
to n> = (3n/2) — ^tan^ 1 {U a b / . For a bounded unperturbed DOS, numerical DMFT calculations give results 
closely paralleling those described above [2(|. At a critical coupling, g < g c = (U a b/2W) c , this line of first-order 
transitions ends in a second order quantum critical end point (QCEP) at T = [2l| . 

Away from A = or t a b = y/t aa tbb, the terms, AJ2i{fLfib + h.c) (first case) or {t aa t bb - t 2 ab ) J2 <l , j> {fLfjb + h.c) 
(second case) change the local U(l) invariance of Hfkm to a global U(l) invariance: only the total fermion number 
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FIG. 1: Local spectral function of the U a b < periodic Anderson model within DMFT in the symmetric (black, dashed) and 
asymmetric (red, dotted) limits for the a (lower panel) and b Fermions (upper panel). Upon varying e' a , the charge Kondo 
insulator (e' a = 0) melts into a charge-Kondo-FL state for e' a 7^ 0. 



commutes with H, [nf a +n,f b ,H] =0. This global symmetry can indeed be (and is) spontaneously broken, and, away 
from "half-filling" , gives a color sing let superfiuid (CSS) at low T, characterized by A Qb = £\ e ab (f}jj b ) > below 
Tsc = fiUab/tfAjTi). This is readily checked by studying H within the generalized Hartree-Fock random phase 
approximation |6|, which is valid now, since there is no local constraint for any A / 0. The enhancement of the 
superconducting T c in doped PbTe and SnTe 0] is now rationalized 11 1 as a manifestation of the proximity to the 
resonant charge fluctuation regime driven by strong quantum fluctuation between the two degenerate valence states, 
11 f b =0,1, envisaged in the original prediction [Io| . for A ~ as above. The respective transition temperatures can 
be computed from DMFT over the full range of \U a b\/ (t aa , tbb, t a b), and yields a smooth BCS-BEC crossover with 
increasing \U a b\- 

This implies that the charge Kondo screening can never produce a CKL ground state. However, above the ordering 
scales, relevance/irrelevance of the /„ — /& hybridization determines whether the IBM or a correlated FL metal will 
result. Consider the symmetry unbroken phase for A > 0. The term A X)-t(/i L/»* + ^- c ) hybridizes the "infinitely 
heavy" (for A = above) fb fermion with the / a -band. This is precisely the U a b < periodic Anderson model (PAM) 
with a local hybridization between the "light" (a) and "heavy" (b) fermions. To access the physics in this regime, we 
solve this PAM within DMFT. 

We choose a Gaussian unperturbed DOS for the a-fermions with W/H — 12kHz and \U a b\/H = 20kHz as represen- 
tative values. To solve the impurity problem of DMFT, we use the iterated perturbation theory (IPT) as an impurity 
solver. Prior work carried out for the U > periodic Anderson model with additional f — d interaction [23j has 
demonstrated its accuracy vis-a-vis the more "exact" QMC solver [22| : the IPT spectral functions closely match those 
extracted from QMC at considerably less numerical cost. Additionally, self-energies can be readily extracted from 
IPT in contrast to the difficulties encountered with QMC in this context. In Fig. [TJ we show the results obtained 
within DMFT for the PAM for two cases: (i) e' a = 0, and (ii) e' a ^ 0. The first case has particle-hole symmetry, while 
the second does not. For a half-filled a-band, in case (i), we find a hybridization-induced Kondo insulator, where the 
hybridization gap opens up after local dynamical correlations have generated a charge-pseudospin quenched analogue 
of the Kondo screened state. The appearance of the high energy Hubbard bands is an explicit manifestation of rel- 
evance of local dynamical correlations. At lower energies, hybridization opens up a (renormalized) band insulating 
gap. Interestingly, this is a new type of insulating state, which we dub a Charge Kondo Insulator (CKI). Obviously, 
since G a b(u>) = (f a ; fl) is finite, this CKI state supports a finite excitonic polarization in the sense of Portengen et 



al. |24|. 



In case (ii), lack of particle-hole symmetry, i.e e' a 7^ 0, effectively leads to partially filled a, b bands, inducing a 
correlated metallic state within DMFT. In this case, the T = spectral function of the b ("heavy") fermion shows 
a sharp quasiparticle resonance near u) = 0, and examination of the self-energy shows that the metallic state is a 
renormalized Fermi liquid (FL). "Related signatures of the CKI to CKL transition are also clearly visible in the 
"conduction" (a) band DOS." Thus, DMFT yields a correlated FL metal due to dynamical quenching of the charge- 
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pseudospin degree of freedom by itinerance of the 6-fermion (due to "hybridization", A, in Eq. (JTJ) ) . Emergence of 
this charge-Kondo-FL (CKL) within DMFT for the U a b < PAM is understood as follows. In the impurity version 
of the lattice problem, we have an X-Ray edge problem with recoil (Eq. @) for A ^ 0. Within DMFT, the infra-red 
singularities described above for A = are immediately cut off below a low energy scale associated with this recoil. 
Expressing the low-energy part of H in terms of collective, bosonic variables pBj . this coherence scale, below which 
correlated FL behavior sets in, is estimated to be ksTcoh — Z ~ exp(£/^ b £n(7)/(l — 7 2 )) with 7 = m,f a /mfb, the mass 
ratios of the f a ,b fermions. In our DMFT calculation, we find Z ~ 0.3 from the real part of the self-energy (not shown). 
Hence, with % = 0.3£ a , we obtain T co h — 0.071>. With a slightly smaller value of U a b/W, T co h can be increased to 
O.lTp, which would lie in the lower range of temperatures observable in current experimental realizations 0. Above 
T co h, the irrelevance of the "hybridization" (A) implies that the system crosses over smoothly to the IBM phase found 
for A = 0. Away from t a b = \Jt aa tbb, the /^-fermions once again acquire a finite mass, giving correlated FL behavior 
below T co h, or broken symmetry ground states. 

Only if T co h > T$c.cdw will the CKL state be observable in a small window in parameter space. The above 
analysis suggests a way to observe the CKL state experimentally in optical lattices: 

(a) force a mixture of two atomic states via a Feshbach resonance (creating attractive interaction) , 

(b) apply sequential, time-synchronized laser excitation to simulate the "hybridization" , and 

(c) "tune" A to relevance in the RG sense, by manipulating the one-body potential [lit . 

The change in the Fermi surface (FS) as the system is tuned across t a b = \fUuJb~b, A = could be detected by 
time-of-flight measurements 26] for T > Tsc,cdw ■ At t a b, A = 0, a single smeared FS, corresponding to itinerant f a 
quasiparticles, should be seen, while, for t a b > y/t aa tbb and/or A^O, two well-defined FS sheets (k-space eigenvalues 
of the free part of H in Eq. |T]), corresponding to coherent quasiparticle propagation, should be observable. Since 
the self-energy is strictly local in DMFT, the FS will be unaffected by interactions, which then have the sole effect 
of critically (or overcritically) damping the FS. Given that the CKL state involves local, dynamical renormalization, 
designing systems where the actual, measured FS agrees well with that computed without interactions will show 
improved chances of experimentally observing the CKL state. 

Recently, the quasiparticle dispersion and the spectral function of strongly interacting fermions exhibiting a 
BCS-BEC crossover in an optical trap has been experimentally deduced by Stewart et al. 27]. We suggest 
that similar experiments carried out below T co h — 0.07Tp should exhibit a relatively sharp low-energy peak in 
Ak(ui) = (— l/7r)Im[l/(cj — €ka — £a(w))] _1 : this would constitute a characteristic signature of the CKL state. The 
CKI state in the symmetric limit above could also be probed using this technique: one should observe two (a and b 
bands) dispersive features below and above Ep, with a finite hybridization gap at Ep. At higher T currently accessi- 
ble, we predict an incoherent pseudogap structure in Ak(u>) at low energy. Achieving lower T may be feasible in the 
near future: theoretical suggestions to do this by efficient shaping of the optical lattice profile is predicted to lower 
the entropy per particle by a factor of 10 [281 ]. 

Our study also permits drawing an analogy with the phase diagram of nuclear matter. In quantum chromodynamics 
(QCD), where three quarks form a trion bound state, a baryon in the N — 3 (color) theory, one also has, as a function 
of chemical potential, "color superfluid" (CSS) and high-T, quark matter phases. In the cold-atom setting, in a 
pioneering work, such a trion gas phase has theoretically been shown to be separated from a color superfluid phase by 
a second order transition [29I ]. Let us now study the possibility of realizing analogous physics in our case. Our model 
corresponds to the "N — 2" version of QCD, where the CSS (p — p order) phase is separated from the excitonic (p— h 
order) and incoherent IBM phases by a line of first order transitions, ending in a second order end-point, generically at 
finite T. Interestingly, the IBM phase has co-existing single-particle ("quark") and pair states at low-energy; it is an 
incoherent state where "quarks" and "pairs" continuously transmute into each other, their equilibrium fraction being 
determined by the ratio U a b/t,A and T. The confinement-deconfinement transition, familiar in the QCD context, 
can be illuminated in our case by employing bosonization of the "impurity model" of DMFT. At A = 0, the impurity 
version of Hp km is readily bosonized [301 ] into a collection of radial bosonic models about the "impurity" . This reads 
Lq = vp J[n 2 (r) + (d r (f>(r)) 2 ]dr, where vp = 2ta is the carrier fermi velocity. For small <5 = (t a b — ^t aa tbb) ^ 0(< 0), 
the f a — ft hybridization is small enough that the two-particle hopping terms (second order in S) are more relevant. 
These terms are H res ~ — (5 2 /U a b)J2<ij>fiafibfj b fj a > which can be decomposed in two ways in a Hartree-Fock 
decoupling (exact in D = 00, and in the corresponding impurity model, this corresponds to two, p — h and p — p 
pairing fields), generating either excitonic (particle-hole) order (in the imbedded impurity problem of DMFT) with 
order parameter A p h = J2i(fLfja + fibfjb) > or CSS (particle-particle) order, i.e, J2i t<Lh (fLfib) > 0- Bosonization 
then yields a dual cosine quantum sine Gordon (QsG) model: 
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L = L - gi / cos(l3i4>(r))dr - g 2 / cos((3 2 0(r))dr (3) 



with Pi — \/8ttK and 02 = W&njK. K = [1 — (t//27rt>F)] _1 in the weak coupling limit for our negative— U model 
exactly at the FK point in the impurity model. A renormalization group analysis of L shows two different phases, 
depending upon the sign of Sg — (gi — g 2 ) [3lj ]: 

(i) when gi > g 2 , the (j)(r) field acquires a finite expectation value, giving p — h excitonic order. In QCD lore, this 
translates into (ipaipa) > 0, with a = a, b, i.e, to a chiral meson phase, and, 

(ii) when gi > gi, the 0{r) field acquires a finite expectation value, giving CSS order. This corresponds to 
(V'aV'b) > 0> i- e i to a color-singlet superfluid in QCD. 

(Hi) Interestingly, these phases are separated by a critical curve gi = g% or K = 1, where the QsG model becomes 
self-dual This is exactly soluble, and the correlation functions are characterized by continuously varying, iGdependcnt 
exponents. 

The above meson and CSS phases arise from the IBM (bosonic gaussian model) via Kosterlitz-Thouless transitions 
(driven by relevance of either of the cosine terms for K < 1 or K > 1), and are separated from each other by a line 
of first-order transitions, ending in a second-order endpoint at finite T = T* . For T > T* , both (tptp), (tptp) = 0, 
corresponding to "quark matter" in QCD [32j]. Of course, this analogy is only suggestive: in our case the a, b particles 
have very different masses, in contrast with the small mass difference between u,d quarks in QCD. Nevertheless, the 
analogy is interesting, and generalization to N = 3 is possible. We will extend this qualitative analogy in more depth 
in future. 

In conclusion, from a DMFT calculation, we propose that a combination of wilfully tunable Feshbach resonances 
and sequentially applied laser pulses in optical lattices to tune U a b/t a b, A through a desired range, could facilitate 
the observation of the elusive CKL state. Time-of-flight measurements can, by probing the "Fermi surface" , provide 
valuable guidance on the experimental conditions favoring observation of the CKL state. Our analysis also shows 



interesting analogies with the "phase diagram" of nuclear matter. Extensions to asymmetric interaction models 14j . 
allowing for realization of more exotic ordered phases and anisotropic CKL states within the DMFT approach, is also 
foreseen. 
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